Magnetic blockade mechanism for quantum nucleation of superconducting 
vort ex- ant i vortex pairs in zero external magnetic field 
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We propose a magnetic dual of the Coulomb blockade effect for quantum nucleation of flux vortex 
pairs in high-Tc superconducting (HTS) films and grain boundaries in zero applied field. The 
magnetic blockade instability occurs at 9 = ir, where 8 is the "vacuum" or theta angle. The 9 
term has recently been discussed in the context of several other systems, including charge and spin 
density waves, topological insulators, the quantum Hall effect, and spontaneous CP violation. Our 
model predicts a sharp pair creation threshold current a.t 9 = n, analogous to the Coulomb blockade 
voltage of a tunnel junction, and explains the observed thickness dependence of critical currents 
in HTS coated conductors. We use the Schrodinger equation to compute the evolving macrostate 
amplitudes, coupled by a generalized tunneling matrix element. The simulations yield excellent 
quantitative agreement with measured voltage-current characteristics of bi-crystal and other HTS 
grain boundary junctions. The model also predicts non-sinusoidal behavior in the voltage oscillations 
resulting from time-correlated vortex tunneling. 

PACS numbers: 03.75.Lm, 74.50.-(-r, 74.25.Wx, 74.78.-w, 74.72.-h, 74.70.Xa 



I. INTRODUCTION 

The last several decades have seen a crumbling of the 
barrier between the quantum and macroscopic realms. 
High-Tc superconductivity provides an exquisite exam- 
ple of long range quantum coherence at temperatures 
above 77 K. Cooperative quantum tunneling has also 
emerged as an important class of phenomena, whose 
manifestations include coherent Josephson tunneling^ 
and tunneling of vortices, which can occur more 
readily in layered high-Tc superconductors due to their 
pancake-like structures <^ 

Vortices are known to nucleate in supcrfluid helium 
above a critical velocity,^ underscoring the point that 
magnetic flux is not crucial to vortex nucleation. In fact, 
as will be discussed, the magnetostatic energy contained 
within a superconducting flux vortex actually impedes 
vortex nucleation for applied currents below a critical 
value. Both Abrikosov and Josephson vortices, the lat- 
ter being especially important in high-Tc superconduct- 
ing (HTS) grain boundaries, play important roles in lim- 
iting the critical currents of HTS films and tapes. A 
Josephson vortex is a 27r flux soliton in a wide Josephson 
junction (wide compared to the Josephson penetration 
length Aj), whose energy includes the extra Josephson 
coupling energy in the region where the phase difference 4> 
across the junction advances by 2tt from one minimum 
to the next. 

Though spatially extended, Josephson vortices are ex- 
tremely light. Using the expression provided by Gros- 
feld and Stern, ^ the effective mass of a Josephson vortex 
spanning a l-fim thick film is estimated to be ^ lO^^me. 
The mass of a single pancake Josephson vortex (Joseph- 
son pancake^) within the /c*'* superconducting cuprate 
or pnictide layer, whose phase difference (j)k advances by 
27r, is several orders of magnitude smaller still. An ap- 



plied current in zero external magnetic field can induce a 
vortex- antivortex pair to nucleate in close proximity by 
quantum tunneling^ and subsequantly expand outward, 
as illustrated in Fig.[IJa). 

Conceptually, this process is analogous to electron-hole 
pair creation via Landau-Zener tunneling through a tilted 
bandgap in a semiconductor or Schwinger pair produc- 
tion of electron-positron pairs out of the vacuum.^'' In 
this case, the applied current / acts as the driving force. 
Neglecting the magnetostatic energy and using the ex- 
pression derived by Cohen and McGady^° for the pair 
production rate, one finds the voltage across a wide grain 
boundary J J to be proportional to I exp[-/o //] . Although 
highly nonlinear, this lacks a sharply defined threshold 
current below which the voltage is identically zero. How- 
ever, recent experiments probing HTS bi-crystal grain 
boundaries with femtovolt precisionii suggest the exis- 
tence of a sharp critical current below which no Joseph- 
son vortices nucleate, either in the middle or at the edges 
of the grain boundary. 



II. MAGNETIC BLOCKADE THRESHOLD 
CURRENT FOR VORTEX NUCLEATION 

The Coulomb blockade effect is well known for sin- 
gle electron tunneling,!^ and has also been proposed for 
charge soliton pair creation in the massive Schwinger 
modelfi^ essentially a quantum sine-Gordon model that 
includes electrostatic effects, and in density wavesJ^"— 
Here we propose a magnetic blockade threshold current 
for Abrikosov or Josephson vortex pair creation^ in HTS 
films and grain boundaries.. -J This threshold, for Joseph- 
son vortex nucleation in a grain boundary, can be much 
smaller than the Ambegaokar-Baratoff critical current 
even when d-wave pairing symmetry is taken in to ac- 
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count .J^ A magnetic "Weber" blockade effect has also re- 
cently been proposed to interpret magneto-conductance 
oscillations of thin, narrow superconducting strips.— 

Figure [IJa) illustrates nucleation, in zero applied mag- 
netic field, of a pair of Josephson vortices in a wide HTS 
grain boundary. The vortex pair generates a magnetic 
field (yellow shading): B* — /3$o/A^ that links the vor- 
tex and the anti- vortex. Here A = vXjXl or A = Al 
for a Josephson or Abrikosov vortex pair, respectively, 
$0 = h/2e is the flux quantum, Al is the London pen- 
etration length, and /3 ^ I/tt is a geometrical factor. A 
current density J in a film of thickness d creates a mag- 
netic field, Bj « /io J(i/2, which can partially cancel the 
internal field linking the vortex and antivortex and thus 
overcome the magnetic barrier (Coulomb barrier analog) , 
provided its value is high enough. 

At low applied currents, the difference in magneto- 
static energy densities with and without the vortex pair, 
(Sj±B*)V2Ato-Sj/2Mo, is positive when \Bj\ < B*/2. 
However, when 9 = 2ttBj/B* > tt, the formerly lowest 
energy state becomes a metastable state or "false vac- 
uum," as illustrated in Fig. [T](b). The phases 0^ tunnel 
coherently into the lower energy well, creating a bubble 
of "true vacuum" bounded by flux solitons (Josephson 
vortices) when 9 exceeds tt. Setting 6 = tt thus yields 
the pair creation threshold current per unit width: 

j° = B*/fio = /3$o/(moA2). (1) 



a c 




FIG. 1. (a). Nucleation of a Josephson vortex- antivortex pair 
(green and red) in an HTS grain boundary. The yellow shad- 
ing indicates the magnetic flux linking the vortices, (b). u{(j}) 
when 9 = 2nE/E* > tv as the phases <l)k{^) tunnel coherently 
into the lower well. (c). Potential energy vs. 9 for (p ~ 27rn. 
(d). Modified resistively shunted junction model, in which 
time-correlated tunneling of Josephson vortices and antivor- 
tices is represented by analogy to time-correlated SET.— 



which yield the vortex pair creation critical current: 



For the case of Josephson vortex nucleation, where 
A^ = AlAj, this yields a scaling with classical critical cur- 
rent density Jq of the form, j^^ cc Jq^^, since Aj oc Jq 
for both planar and sandwich-type junctionsi^ For cur- 
rents below the vortex nucleation threshold, 9 can be 
written in terms of the supercurrent per unit width js 
as = TT {js/jpc)- If the vortex and antivortex are sepa- 
rated by a distance then the pair creation critical cur- 
rent between vortices can be written as jp^,^ ^ ^o/2L, 
where L ~ /xoA^/2/3f is the inductance of the fiux toroid 
coupling the vortices. This is essentially the dual of the 
Coulomb blockade voltage Vc — e/2C for a small capaci- 
tance tunnel junction. Nucleation of a vortex near one 
edge and antivortex near the other can be treated simi- 
larly using the image vortex concept.™ 

HTS coated conductors generally have small-angle 
grain boundaries^ and their critical currents vs. thick- 
ness d often plateau, so that Jc decreases with d.^^ This 
is readily explained by the vortex pair creation current 
density, Jp^ = (3^o/{fioX'^d) when d > X, resulting from 
the magnetic blockade effect. For a strip of width w, 
the pair creation current becomes: 1°^ = f3^ow/{fio^'^), 
which is independent of d until it becomes large enough 
for vortex rings to nucleate. When d << A, the size of 
the vortices increases since the effective 2-D penetration 
length A scales inversely with d:— A = A^/d. Since A = A 
when d >> A, we use the approximation: A sa A -I- /d, 



Figure [2] shows a favorable comparison of Eq. ([2]) 
with measured critical currents vs. thickness of HTS 
coated conductors. The introduction of multiple insulat- 
ing Ce02 spacer layers into HTS coated conductors has 
been found to increase critical current. Our model pro- 
vides a straightforward interpretation - sufficiently thick 
spacer layers decouple vortex pair nucleation events on 
adjacent HTS layers, enabling Ic to increase with the 
number of layers. 



1000 




Thickness (|jm) 

FIG. 2. Theory vs. experimenlsS^ for critical currents of 
HTS coated conductors vs. thickness. (See text: we use A= 
380 nm, J°^= 1030 A and 398 A to fit 2008 and 2005 data, 
respectively.) 
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III. TIME CORRELATED VORTEX 
TUNNELING: THE SCHODINGER EQUATION 
AS AN EMERGENT CLASSICAL DESCRIPTION 

The "vacuum" or theta angle 9, discussed extensively 
in the quantum field theory literature (e.g. see Ref.-iS, and 
citing papers) is related in our model to total displace- 
ment flux $ by: 9 = 27r($/$o), where $ = X'^Bj//3 + 
n^Q. The potential energy of the fc*'' layer for a cuprate 
or pnictide superconductor can then be written, similar 
to the massive Schwinger model, ^'^ as: 

u[(l)k] = Uj [1 - cos4ik{x)] + {B - (l)k{x)f , (3) 

where the first term is the Josephson coupling energy and 
the quadratic term is the magnetostatic contribution. 

Figure [ijc) shows plots of u vs. 9 when the energy 
is minimized for 27rn (dropping the subscript k) 
when um << uq. Each u vs. <j) branch has a parabolic 
form, oc (0 — 2?™)'^ cx (<!> — n$o)^/2i', analogous 
to the charging energy parabolas, oc (Q — ne)^/2C, in 
time-correlated SET^^ and time-correlated soliton tun- 
nelingji^ii^ The phases 0fe tunnel coherently into the next 
well as each parabola, or branch, in Fig. [Ijc) crosses 
the next at = tt and the other instability points 
9 = 27r(n + l/2). Regardless of the detailed shape of 
the Josephson coupling energy (which may or may not 
be sinusoidal and may even include disorder), the behav- 
ior that emerges is clearly nonsinusoidal in nature, which 
affects the Shapiro steps in the I—V curve in the presence 
of a microwave signal. 

The approach proposed here is based on the hypoth- 
esis that the amplitudes tpn and Tpn+i for the system to 
be on branches n and n -I- 1 , respectively, are coupled via 
coherent tunneling of Josephson pancakes, represented as 
microscopic quantum solitons.^'^'^- This picture is moti- 
vated by Feynman's intuitive derivational of the dc and 
ac Josephson effects. Advancing (f>k{x) by 27r in a fi- 
nite region is equivalent to creating a pair of microscopic 
solitons - essentially pancake Josephson vortices^ with 
extremely small effective massesil Similarly, Abrikosov 
pancake vortices^ are pointed out to have small effective 
masseSi^ In this picture, many coupled pancake vortices 
behave as a quantum fluid, within which their coherent 
tunneling is viewed as a secondary Josephson effect. 

We use a modified version of the resistively shunted 
junction (RSJ) model [Fig. [lid)], in which the analogy to 
time-correlated single electron tunneling is used to simu- 
late time-correlated vortex tunneling and obtain voltage- 
current characteristics of an HTS grain boundary junc- 
tion. The voltage across the junction is: V — d^/dt — 
{^o/2Tr)d9/dt, yielding a normal current through the 
shunt resistance: /„ = V/R — {^o/2TrR)d9/dt. For a 
grain boundary junction with total current /, we define 
LJ = (27ri?/$o)/- 

In a wide junction, advancing the phase by (0) (the av- 
erage phase for all layers) near the middle of the junction 
creates a Josephson kink-antikink pair whose circulating 



currents either reduce or reinforce the existing supercur- 
rent between kinks. This yields a net supercurrent given 
by: Is = (/*/27r) [9 - (0)]. Using /„ = / - Is and intro- 
ducing the time constant r — L/R, yields the following 
equation for the time evolution of 9: 

df) 1 

^ = .^-[9^m- (4) 

The expectation value (0) is computed by solving the 
time-dependent Schrodinger equation: 

= UaipOA + JVV'i.o , (5) 

viewed as an emergent classical equation following Feyn- 
maufSI to compute the original and emerging macrostate 
amplitudes ipo{t) and ipiit)- These are coupled via a 
tunneling matrix element with a Zener-like force de- 
pendence, which transfer the phase cjjk from one well to 
the next [Fig. [1Kb)]. 

Our model represents the amplitudes ipo^i for the sys- 
tem to be on branches and 1 [Fig. [IJc), more generally 
n and n -|- 1] by: 

"00,1 = ^po,i expi(5o,i . (6) 

For a layered superconductor, po.i — Nq^i/N is the frac- 
tion of superconducting layers on the respective branch. 
In this picture ?/'o,i are not this superconducting order 
parameters on opposite sides of the junction but, rather, 
represent the amplitudes for the phases (t>k to be in ei- 
ther of the two wells shown in Fig.[TJb). The phases Jq.i 
are the Berry phases for the vortices, which would, for 
example, lead to an Aharonov-Casher effect for vortices 
traveling along two branches of a ring-shaped Josephson 
junction surrounding an island charge.— 

Advancing 4>k{x) by 2tt within a given region, taking 
(j)k from one branch to the next, is equivalent to creating 
a pair of microscopic 27r-solitons (Josephson vortices). 
The driving force is the energy difference per unit length 
between potential minima at (/) ^ 2?™ and (f> ^ 27r(n+ 1). 
When um/uq << 1, this force is: F = 47rujv/6'„', where 
9J = 9-2n{n+^). 

Following Bardeen's procedure^^ii^ the matrix ele- 
ment for Zener tunneling through the soliton energy gap 
is estimated as: 

T^(F) = -4^^Aexp[-Fo/F] , (7) 

where ^ A^p/hvo + A^^, Am is a mean free path 
length, Aip is the microscopic soliton energy for a pan- 
cake Josephson vortex, vo is the phase (Swihart) veloc- 
ity,— and Fg ~ A'^/hvo- This expression for is similar 
to the rate of Zener tunneling or Schwinger pair produc- 
tion in l-D.-ifl Since any negative energy difference within 
the "bubble" is balanced by the positive flux soliton pair 
energy, the matrix element couples states of equal energy, 
Uq = Ui = U. Thus, defining t/^o,! = xo,i(i) exp[~*C^V^]j 
the Schrodinger equation [Eq. ([5|)] reduces to: 

ihdxos/dt = T^xi,o- (8) 
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Below the pair creation threshold, we can write: 6 = 
21:1/1* where /* = 2/3$o^/(AioA^) = {^o/L){t/w) within 
a grain boundary region of width £, which may be 
slightly smaller than w due to the finite width of the 
vortices. Due to multiple vortex pair nucleation events 
in a wide junction, H/w may not decrease monotonically 
with w/Xj. The variables are put into dimensionless 
form: t' = t/r, f = a;r/27r, proportional to total cur- 
rent I, q = 9/2tt, where qo = 9o/2n = Fo/(47ruAf), 
q„' = 6'„'/27r = q — n— 1/2, and is = Ig/I* = q — p—n is 
the normalized supercurrent, and the normalized voltage 
i^- " = ^ = / ^ *«■ Finally, setting xo(0 = co{t) and 
Xi(i) — ici(t) in Eq. (jH)), taking cq and ci to be real, 
yields the following coupled difi^erential equations: 

-177 ^ [7'7n'exp(-go/<?r/)]co , 
dco 

= -[7'7«'exp(-(7o/9n')]ci , 

for > and where 7 = 327r^UMA. These are inte- 
grated numerically with initial values cq = 1 and ci = 0. 
The phase expectation value is: {4>) = 27r[n-f p], where 
p = I ci p , and the transition is considered complete when 
p exceeds a cutoff close to one. 

For each value of / — Ijl*-, a time average over 
several complete cycles is performed to compute v — 
[l I w){y 1 1* B}) . A similar approach has recently been 
used by the authors to compute current-voltage charac- 
teristics due to soliton pair production in density waves.— 
As in the density wave case, each transition from one 
branch to the next takes place over a relatively long time 
interval, suggesting that many pancake vortices "flow" 
through the barrier like a quantum fluid rather than tun- 
neling abruptly as a single massive flux line. 

Figure [3] shows resulting theoretical normalized volt- 
age vs. current plots for several values of 7 for a fixed 
value of go =3 and Ijw — 1. The V-1 plots exhibit 
piecewise linear behavior, which fits neither the classical 
RS J model nor a straightforward thermally activated flux 
flow model. Similar piecewise quasi-linear behavior is 
seen quite frequently in high- Tc superconducting YBCO 
grain boundary junctions. The main effect of increasing 
7 is to increase the slope of the Y-1 curve, saturating at 
a slope of one for the normalized plot. 

Figure m shows similar plots for a fixed value of 7 = 20 
and several values of (Jq. Here, the main effect of in- 
creasing which corresponds to increasing the energy 
required to nucleate each vortex-antivortex pair, is to in- 
crease the degree of rounding in the V-1 plot. Moreover, 
note that the simulated "measured" critical current 1^ 
increases to well above the nominal magnetic blockade 
critical current /pc ^ /*/2 as go becomes significantly 
greater than one. 

Figure [SJa) shows comparisons between theory and 
the measured V-1 characteristics of a YBCO^ grain 
boundary junction at several temperatures. The sim- 
ulated 86 K plot in Fig. [SJa) (top) was obtained us- 
ing the classical RSJ model in the overdamped limit: 




2 4 6 8 
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FIG. 3. (Color online.) Normalized theoretical voltage- 
current plots for a fixed value of 50 = 3 and several values 
of 7. 




FIG. 4. (Color online.) Normalized simulated voltage-current 
plots for a fixed value of 7 = 20 and several values of go. 



V/{lcRn) = \/ {I /IcY — 1, without invoking thermal ac- 
tivation. The fact that the 86 K data fits the classical 
RSJ model almost perfectly suggests that the effective 
Josephson penetration length is comparable to or longer 
than the junction width (short junction limit) due to the 
small Josephson coupling energy at this temperature. 

The theoretical plots for the remaining temperatures in 
Fig.IHKa), as well as for Figs.[5jb) and[5Kc), were obtained 
from the Schrodinger equation as discussed above, using 
the parameters shown in Table HI The increase in q^ going 
from 82.5 K down to 70.0 K is consistent with increasing 
Josephson coupling energy, leading to a larger _Fo E^nd 
smaller Aj, as the temperature decreases. The increased 
rounding of the V-1 curves with decreasing T is also con- 
sistent with the system going from the short- to the long- 
junction limit as Aj decreases. Moreover, the fact that 
the V-1 curves are more rounded for low than for high 
temperatures provides powerful evidence that the round- 
ing of the V-1 curves is not due to thermal activation. 
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The extremely light masses of Josephson vortices coupled 
with quantmxi fluidic properties, discussed above, enable 
quantum effects to dominate over the entire temperature 
range. Further evidence that even Abrikosov vortices ex- 
hibit quantum, rather than thermally activated, behavior 
over a wide temperature range is provided by magnetic 
relaxation rates of trapped flux in YBCO, which go down 
with increasing temperature^ from 4 K to about 86 K. 




FIG. 5. (Color online.) (a). Experiment vs. theory for a 
YBCO bicrystal junction.— Top curves: V-I curve at 86 K 
(dotted line) vs. classical RSJ model (solid line). Remain- 
ing curves: quantum nucleation model (solid lines) vs. mea- 
sured V-I curves (dotted lines) at (top to bottom): 82.5 K, 
77.2 K, 75.0 K, and 70.0 K. (b). Comparison between quan- 
tum nucleation model (solid line) vs. V-I characteristic 
(diamonds) of an iron pnictide superconductor bicrystal at 
4.2 K.— (c). Model simulation vs. experiment for V-I curve 
of grain boundary junction in a thallium-based cuprate su- 
perconducting film at 77 K.— 



TABLE I. Parameters used to generate the simulated V-I 
curves in Fig. [S] 



Figure Plot 


ic/r 


w/l 


7 


qo 


Fig.EJa) - 82.5 K 


0.55 


1.15 


50 


0.5 


Fig. IHa) - 77.2 K 


0.72 


1.09 


45 


1.7 


Fig. [Ha) - 75.0 K 


1.25 


1.16 


40 


6.5 


Fig. IHa) - 70.0 K 


1.80 


1.16 


30 


12.0 


Fig.IHb) 


0.71 


1.04 


60 


1.9 


Fig. [He) 


1.00 


1.00 


12 


6.0 



and generalized tunneling matrix element T^, coupled 
with the analogy to time-correlated single electron tun- 
neling suggested in Fig. [TJ provides a simple, yet pow- 
erful approach to modeling the dynamics of vortex tun- 
neling. The top plot in Fig. [HI shows predicted voltage 
oscillations, exhibiting highly non-sinusoidal behavior, 
for a uniform junction with the parameters 7 = 40 and 
go = 6.5, for an applied current / = ///* = 2.5. 
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FIG. 6. Simulated normalized voltage oscillations vs. nor- 
malized time for (top-to-bottom, offset for clarity): a uniform 
junction with 7 = 40 and go = 6.5; and nonuniform junc- 
tions, represented as 100 junctions in parallel using the same 
average 7 and go, but (see text) with 5 = 0.1, k — 0.05, 
S = 0.1, k = 0.1, and S = 0.5, k = 0.1. 



Figure EJb) shows excellent agreement between our 
quantum model and the V-I curve of an iron pnic- 
tide superconducting bicrystal, '^^ consisting of coupled 
SrFei.74Coo.26As2 and Bao.23Ko.77Fe2As2 crystals, each 
~300 fim wide, and thus in the long junction limit. Sim- 
ilar piecewise linear behavior to that seen in the figure 
also occurs frequently in cuprate grain boundary junc- 
tions [e.g. Fig. IHa), 77.2 K data and Ref.37.38] ^nd is 
analogous to the piecewise linear /- V curve of an ideal 
Coulomb blockade tunnel junction. Figure ^c) shows a 
plot of our simulation as compared to the measured V- 
I characteristic of a thallium-based cuprate,— exhibiting 
good agreement with the data within experimental error. 
Table U] shows the parameters used for the simulations. 

The use of the time-dependent Schrodinger equation 



In the remaining plots of Fig. [S] we model non- 
uniformities by representing the junction as 100 junc- 
tions in parallel, for which 7 and qo are psuedo-randomly 
varied with a uniform distribution within the ranges 
7(1 ± S) and qo{l ± S), respectively. The phases be- 
tween adjacent domains i are coupled via an addi- 
tional term that adds to the net force for each domain: 
/c [((/) — « + 1) —2 applying periodic bound- 

ary conditions. The remaining plots in Fig.lHlalso assume 
f — 2.5 and the same average values for 7 and go as be- 
fore, but with different values of 6 and k. The middle two 
plots, for which i5 = 0.1 and k = 0.05 (second from top) 
and k — 0.1, show greatly reduced voltage oscillation am- 
plitudes as well an apparent amplitude modulation effect 
consistent with reports of subharmonic Shapiro steps4S 
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The bottom plot {5 = 0.5 and k = 0.1), representing the 
highest degree of disorder for this series, shows osciUation 
amphtudes that are even further reduced. This is con- 
sistent with the idea that Shapiro steps would essentially 
disappear in a highly nonuniform junction. 

IV. CONCLUDING REMARKS 

A key premise of this paper is that the concept 
of coherent Josephson-like tunneling should be gener- 
alized to include other phenomena with quantum flu- 
idic properties. Rather than treating the system as a 
massive object and computing the total Euclidean ac- 
tion or "bounce," a generalized tunneling Hamiltonian 
matrix element is introduced to connect the original 
and emerging macrostates. The "disconnectivity" be- 
tween macrostates, introduced by Leggett,^^ is extremely 
small as in Josephson tunneling. This notion is further 
supported by examination of Schrieffer's real-space de- 
scription of the BCS ground state:*^ V'o = A(po{ri — 
r2)xi2---'y5o(rAr_i - rN)xN~i,N, where ipo is the relative 
coordinate wavefunction of a pair, x is the corresponding 



spin function, and A is the antisymmetrization opera- 
tor. As ipo evolves into a state, ipvpj containing a vortex- 
antivortex pair, each pair state individually evolves into 
a deformed state ifyp topologically connected to tpo- 

Beyond improving our understanding of superconduct- 
ing flux vortex nucleation and dynamics, the concepts 
proposed here could have scientific impact in other ar- 
eas, such as 6* = TT instabilities in spontaneous CP viola- 
tionj^ charge and spin density wavesj^ topological insu- 
latorSj^!^ and the quantum Hall effect4S Finally, better 
understanding of the quantum behavior of flux solitons 
could potentially lead to topologically robust forms of 
quantum information processing. 
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